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1.  INTRODUCTION  AND  BACKGROUND 


This  report  describes  a  formulation  developed  for  analysis  of  radiation 
from  slot  antennas  and  off-surface  wire  radiators  on  open  or  closed  asymmetric 
bodies,  denoted  here  as  bodies  of  translation  (BOT).  Representative  examples 
of  such  configurations  are  shown  in  Figure  1.  The  case  of  near-fields  in  the 
vicinity  of  such  antennas  is  treated  as  a  subcase.  The  formulation,  based  on 
the  method  of  moments  (MM),  is  a  generalization  of  the  results  given  in 
References  1-3.  The  computer  algorithm  A-STAR  (Algorithm  for  jJurface  of 
Translation  Attached  Radiators)  implements  this  analysis  and  is  described  in 
detail  in  Vol.  II  of  this  report.  The  listing  of  the  codes  is  given  in  Vol. 
III. 

The  >01  technique  has  been  applied  to  a  plethora  of  electromagnetic 

problems  that  are  associated  with  complex  aerospace  systems.  Generally,  the 

classical  partial  differential  equation  formulations  of  such  problems  are 

analytically  and  computationally  intractable.  On  the  other  hand,  the 

equivalent  electric-  or  magnetic-field  integral  equation  (EFIE  or  MFIE) 

formulations  of  these  cases  are  often  amenable  to  solution  by  the  MM 

approach.  In  these  integral  formulations,  the  radiating  or  scattering 

structures  and  surfaces  can  be  represented  by  wires,  surface  patches,  or  wire 

grids. The  computational  requirements  of  the  method  to  date  have 

limited  application  of  these  methods  to  bodies  with  surface  areas  on  the  order 

2  2 

of  a  few  square  wavelengths  (A  ).  Larger  surfaces  (of  ~  45A  )  can  be  analyzed 
via  the  MM  technique  if  the  vehicle  body  has  some  degree  of  symmetry,  such  as 
in  the  case  of  bodies  of  revolution  (BOR).  The  present  analytical  technique, 
denoted  as  MM/BOT,  combines  many  of  the  computationally  cost-effective 
features  of  the  MM/BOR  analysis-*  with  some  of  the  shape  flexibility  of  the 
wire-grid  approach.^  This  combination  permits  treatment  of  difficult  boundary 
conditions  associated  with  realistic  radiating  and  scattering  geometries,  such 
as  parts  of  wing  sections  or  noncircular  aircraft  fuselages  with  attached  or 
embedded  antennas. 


Asymmetric  Surfaces  -  Radiation  Analysis 


Asymmetric  Surfaces  -  Scattering  Analysis 


Figure  1.  Body  of  translation  (BOT)  configurations. 
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2.  SUMMARY  OF  COMPLETED  EFFORT 


The  principal  results  of  the  completed  effort  are  summarized  below. 
Detailed  discussions  of  individual  topics  are  given  in  the  sections  indicated. 

•  The  MM/BOT  formulation  was  developed  for  a  general  BOT  configuration 
to  treat  the  radiation  from  asymmetric  aperture  (slot)  antennas  and 
off-surface  radiators  (Sections  4  and  5).  For  radiating  apertures, 
the  formulation  allows  arbitrary  polarization  and  antenna  excitation 
to  be  specified.  The  results  are  shown  to  be  in  agreement  with  the 
predictions  of  the  MM/BOR  analysis-*  and  with  experimental  measurements 
for  wire  antennas  on  surfaces  of  revolution. 

•  The  MM/BOT  formulation  was  specialized  to  compute  all  electric  and 
magnetic  near-field  components  at  an  arbitrary  point  in  the  vicinity 
of  a  radiating  or  scattering  BOT  surface  (Section  6).  The  coupling 
between  wire  and  aperture  antennas  was  treated  (Section  6). 

•  The  computer  algorithm  A-STAR  was  developed  to  implement  all  aspects 
of  the  MM/BOT  formulation.  The  codes,  written  in  FORTRAN  IV,  are 
modular,  machine  independent,  and  user  oriented.  All  parts  of  the 
codes  were  tested,  validated  (Section  7),  and  installed  on  the  RADC 
computer  system.  A  user/systems  manual  (Volume  II  of  this  report) 
describes  in  detail  the  workings  of  the  algorithm.  A  series  of 
example  problems  is  provided  to  illustrate  the  computational  procedure 
for  the  prospective  user. 

The  foregoing  results  have  special  relevance  to  RADC  in  the  following 
areas : 

•  The  generalized  MM/BOT  formulation  allows  efficient  representation  of 
asymmetric  surfaces,  such  as  wing  sections  and  noncircular  parts  of 
aircraft  fuselages,  for  electromagnetic  analysis. 

•  The  analysis  models  both  aperture  and  off-surface  wire  antennas 
located  on  modern  aircraft.  Antenna  radiation  patterns  can  be 
computed  efficiently  for  all  polarizations. 

•  The  near-field  analysis  allows  the  mapping  of  fields  and  currents  in 
the  neighborhood  of  radiators.  The  coupling  analysis  allows 
engineering  criteria  to  be  developed  for  suitable  antenna  placement  on 
aircraft  surfaces. 


•  The  formulation  demonstrates  hybriding  of  surface,  wire,  and  patch 
modeling  techniques  in  a  consistent  analytical  framework  leading  to 
maximum  computational  flexibility  for  EM  problems.  This  framework  is 
adaptable  to  other  EM  codes  developed  by  RADC  such  as  GEMACS  (General 
Electromagnetic  Model  for  the  Analysis  of  Complex  Systems). 
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3.  SUMMARY  OF  PREVIOUS  WORK 


The  earliest  application  of  the  method  of  moments  (MM)  was  for  the 
analysis  of  thin-wire  radiators  and  scatterers. Recent  investiga- 
tors^*^-^  have  used  the  MM  technique  to  examine  such  radiators  on  or  near 
bodies  of  revolution  (BOR).  In  many  real-world  situations  involving  aircraft 
and  spacecraft,  the  principal  physical  features  of  the  antenna  platform 
adjacent  to  the  antenna  feed  point  often  can  be  modeled  as  a  BOR.  Often  the 
electrical  size  of  the  body  to  which  the  antenna  is  attached  Is  comparable  to 
a  few  wavelengths  (X).  In  these  cases,  the  entire  structure  is  coupled 
closely  to  the  antenna  elements,  so  the  whole  configuration  becomes  a 
radiator.  To  analyze  such  configurations  requires  Incorporation  of  complex 
boundary  conditions  in  the  electromagnetic  formulation  used.  The  MM  technique 
is  the  most  appropriate  analytical  tool  for  such  problems. 

To  extend  the  usefulness  of  the  MM  analysis  for  a  broader  class  of 
radiating  geometries,  we  developed  a  MM  formulation  for  antennas  attached  to, 
near  to,  or  embedded  In  bodies  of  translation  (BOT).  Subcases  of  such 
surfaces  are  finite  cylindrical  bodies  of  arbitrary  cross  section,  which  can 
be  open  or  closed.  Previously,  a  MM  formulation  was  used  to  treat  asymmetric 
bodies  using  a  wire  representation  for  the  surface.^  However,  the  computer 
resources  necessary  to  carry  out  the  latter  analysis  are  often  prohibitive  and 
limit  its  general  use.  In  our  formulation,  an  efficient  analytical  treatment 
is  obtained  allowing  substantial  configurational  flexibility  for  both  the 
antennas  and  the  surfaces  on  which  they  are  located.  The  theoretical 
development  parallels  in  part  the  MM/BOR  and  MM/BOR-wire  formulations  in 
References  5  and  3  and  retains  the  modal  expansion  concept  for  the  currents  on 
the  surface. 

The  primary  motivation  for  our  investigation  was  the  efficient  analytical 
treatment  of  antennas  on  or  near  BOTs.  The  results  of  this  work  are  easily 
extended  to  the  scattering  from  such  configurations.  (For  examples  of  the 
MM/BOR-wire  formulation  applied  to  scattering,  see  Reference  17.)  Earlier 
investigators,  such  as  Andreasen,*®  Wallenberg  and  Harrington,^  and  Wilton 
and  Mittra,^  primarily  addressed  the  scattering  from  cylinders  of  arbitrary 
cross  section  and  infinite  length.  The  case  of  finite-length  cylinders  of 
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arbitrary  cross  sections  has  not  been  treated  previously.  Several 
Investigators  examined  the  special  case  of  finite-length,  right-circular 
cylinders  for  various  limiting  cases.  For  example,  Ufimtsev,^1  Kieburtz,^ 
and  Fialkovskii^  developed  solutions  for  thin  cylinders  with  ka  «  1,  while 
Adey^  considered  long  cylinders  with  ka  ■  1,  where  ka  *  2wa/A  and  a  is  the 
cylinder  radius.  Williams^  studied  the  diffraction  from  finite-length, 
hollow  cylinders  where  the  open  ends  did  not  materially  influence  the 
diffracted  waves.  A  complete  study  of  tubular  cylinders  was  made  by  Kao^»^ 
for  arbitrary  length  and  ka,  but  it  was  restricted  to  broadside  illumination. 
His  formulation  resulted  in  a  pair  of  decoupled  integral  equations  that  was 
solved  to  yield  the  axial  and  the  circumferential  currents  on  the  cylinder. 
Recently,  Davis  and  Mittra  examined  the  current  distribution  on  an  open 
cylinder  of  IX  length  and  ka  -  1  using  a  hybrid  formulation  Incorporating  the 
EF1E  and  MFIE  representations  of  Maxwell's  equations.  The  present  Mt/BOT 
formulation  treats  the  foregoing  problems  as  special  subcases. 
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4.  GENERALIZED  MH/BOT  FORMULATION 

In  Che  following  analysis,  as  in  Reference  1,  the  electromagnetic 
radiated  or  scattered  fields  are  expressed  in  terms  of  potentials  from  which  a 
general  integral  equation  for  the  surface  currents  on  the  radiating  or 
scattering  body  can  be  obtained.  For  maximum  generality  and  unity  of 
treatment,  the  EFIE  formulation  is  chosen  over  the  MFIE  formulation  for  the 
currents  on  both  the  off-surface  wire  radiators  and  the  BOT  surface.  The 
generic  configuration  addressed  here  is  shown  in  Figure  2.  The  overall 
(closed)  body  is  composed  of  a  surface  formed  by  an  arbitrary  generating  curve 
translated  along  the  z-axis,  denoted  as  the  BOT.  The  BOT  can  be  terminated  by 
planar  end-caps  with  tangents  normal  to  the  z-axis.  The  assumption  of  a 
closed  surface  and  planar  caps  is  Introduced  to  simplify  the  subsequent 
discussion  and  can  be  removed  if  desired.  The  vicinity  of  the  antenna 
attachment  point  (Figure  2)  is  designated  as  the  junction  region,  consisting  a 
small  annular  disk  on  the  BOT  and  a  wire  attachment  segment  on  the  antenna 
nearest  the  BOT  surface. 


Figure  2  Generic  BOT-wire  configuration. 
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BOT  surface 


Junction 


Figure  2c.  Detail  of  Junction  region  representation. 


4.1  Electric-Field  Integral  Equation 


Imposing  the  usual  boundary  conditions  on  the  perfectly  conducting 
surfaces  in  Figure  2,  the  generalized  EFIE  expressed  in  terms  of  integro- 
differential  operators  L(*)  on  the  BOT,  caps,  wire,  junction,  and  edge  region 
is  given  by 

(E*)  *  -  <ES  )  -  L  (J)  +  L  (J)  +  L  (J)  +  L  (J)  +  L  (J)  ,  (1) 

tan  tan  s  c  w  j  e 

where  (£*)  and  (£s)  denote  the  tangential  components  of  the  impressed 
tan  tan 

and  scattered  electric  fields,  respectively,  on  each  of  the  various  surfaces 
(or  regions)  in  Figure  2. 

The  operator  L  refers  to  the  BOT  surface  S„  and  is  given  by 


Ls(J)  ■  jwu  ffuds-^Vs  ff  (V'vS).  ds  , 


where  the  free-space  Green's  function  is  4  -  l/4xR  exp(-jkR),  J  denotes  the 
current  density  on  Ss,  R  is  the  distance  from  the  source  to  the  field  point, 

Vg  is  the  surface  gradient  on  the  BOT,  u>  is  the  radian  frequency  (sinusoidal 
excitation  is  assumed),  and  v  and  e  are  the  permeability  and  permittivity  of 
the  medium,  respectively.  For  the  caps,  the  operator  Lc(»)  is  identical  to 
Lg(*)  except  the  domain  of  integration  spans  the  caps.  Similarly,  the 
operator  l*e(*)  is  identical  to  Lg(«)  except  the  domain  of  integration  spans 
the  edge  region  on  either  the  BOT  or  the  caps.  Analogously  for  the  wire  part, 


lu(J)  -  Jw  /  J*  ^  ib-  /  si?  • «»  . 


where  Lw(l)  is  the  one-dimensional  operator  on  the  wire  current  density  and 

is  a  unit  vector  along  the  wire.  The  operator  in  the  junction  region 

L.(«)  is  L  (J),  where  the  domain  of  integration  is  restricted  to  the  attach- 
j  w 

ment  segment  Sa;  Lj(*)  when  the  surface  integration  is  confined  to 

the  disk  part  of  the  junction  region  (Figure  2). 

To  solve  Equation  (1)  for  the  unknown  current  density  J  on  the  entire  BOT 
configuration,  explicit  expressions  must  be  obtained  for  all  operators.  We 
use  the  Galerkln  (MM)  technique  to  expand  the  currents  specific  to  each  region 
where  the  curved  BOT  surface  is  subdivided  into  axial  strips,  the  caps  into 
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trapezoidal  patches,  the  wire  radiators  into  connected  straight-wire  segments, 
the  junction  region  into  an  attachment  and  a  disk  part,  and  the  edge  region 
into  patches  on  the  BOT  and  the  caps.  A  set  of  basis  and  testing  functions 
for  the  unknown  currents  in  each  of  these  regions  is  introduced.  The  basis 
sets  for  the  wires  and  the  junctions  are  identical  to  those  given  in  Reference 
3. 


4.2  Basis  Functions  for  the  Currents 

First,  the  expansion  of  the  currents  on  the  BOT  surface  is  considered. 

The  two  orthogonal  components  of  the  surface  currents  on  the  BOT,  i.e., 
j  and  J  ,  are  expanded  in  terms  of  a  modal  expansion  along  z  and  in 
overlapping  triangle  functions  along  the  t-direction.  Specifically,  the 
currents  on  the  i-th  axial  segment  of  the  BOT  are 

K  “  ]£  (In l  5n l  +  rni  ^ni)  <*“P«es/»eter)  .  (2) 

n 


where 


J®“  -  u  f“(t)  v“(z)  ,  a  -  t  or  z 
nl  a  l  n 


(2a) 


and 


vS<*> 


exp(  jnirz/L)  ,  a  ■  t 

exp( jmrz/L)  -  (-l)n  ,  a  -  z 


(2b) 


The  modal  expansion  term  va(z)  spans  all  integer  n  values  from  -N  to  +N. 

n 

l8,  and  1 88  represent  the  orthogonal  current  components  in  the  t-  and  z- 
ni  ni 

directions,  respectively,  corresponding  to  the  i-th  surface  triangle  function 
f®(t),  where 


fj(') 


|t'|,  |t'|  <  1 

,  t'  -  t  -  t.  . 

0  ,  |t'|  >  1  J 


(2c) 


(The  axial  component  of  current  388  is  chosen  to  vanish  at  the  ends  of  the 
BOT.) 
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On  Che  caps,  Che  currenCs  are  expanded  in  Cerms  of  Cwo  componenCs  along 
ut  and  Up,  where  p  is  a  polar  coordinaCe  on  Che  caps  (Figure  3),  yielding  for 
Che  i-Ch  crapezoidal  paCch 


+C  ♦  TCC _CC  .  ♦  TCP,Cp 

JJC  “  “c1*  Ji  +  “p1!  V  » 


where 


J^a  -  f“(c)g“(p)  ,  a  -  C  or  p  (3a 

and  f“(c)  and  g“(p)  are  criangle  funecions  spanning  Che  cap  surfaces.  (In 

X  X 

general,  u  *u  *  0. )  For  maximum  flexlbilicy  in  surface  represenCaClon,  Che 

c  p 

pacchea  are  allowed  Co  be  nonuniform  so  chaC  Che  criangle  funcdons  can  vary 
from  paCch  Co  paCch.  (For  alCernaCe  paCch  represenCaCions,  see  References  29 
and  30.) 


Figure  3.  Detail  of  cap  coordinate  geometry. 


The  basis  functions  for  the  wire  and  the  junction  currents  are  those 
given  in  Reference  3,  and  for  completeness,  they  are  repeated  here.  On 
the  4-th  wire  segment,  the  filamentary  wire  current  can  be  expressed  as 


3"  -  u^I^T^(h)  (amperes)  , 


(4) 


where  u^  is  a  unit  vector  along  the  i-th  segment,  T^(h)  denotes  a  triangle 
function  (Figure  2a),  and  i”  is  the  unknown  wire-current  coefficient 
associated  with  the  4-th  segment.  For  a  wire  segment  nearest  the  attachment 
point,  one-half  of  a  wire  triangle  function  overlaps  the  half-triangle  basis 
term  on  the  attachment  segment.  This  representation  yields  a  piece-wise 
continuous  (p.c.)  current  from  the  junction  to  the  wire. 

The  junction  current  3^,  in  terms  of  basis  functions  associated  with  the 
wire  attachment  segment  Sa  and  the  disk  region  S^,  is  given  as 

♦  1  1  I  ^a  *  p  e  ^a 

JJ(p)  ■  IJ  J  .  (amperes/meter)  ,  (5) 

Pd’  P  e  Sd 

where 


a 

♦  Vh) 

ua  2na 

+  1 

M  —II  - 

Jd 

r  2irr 

-► 

u 

a 

and  T  (h) 

(5a) 

(5b) 


and  ua  and  Ta(h)  are  an  outward-directed  unit  vector  and  a  half-triangle 
function  on  the  attachment  segment,  respectively;  uf  is  an  unit  vector  on  the 
(annular)  disk  surface  away  from  the  wire,  r  is  the  radial  distance  on  the 
disk,  b  is  the  outer  disk  radius,  a  is  the  wire  radius,  and  iJ  is  the  unknown 
junction  current  coefficient. 

The  edge  region  spans  the  rims  of  the  caps  (Sc)  and  the  ends  of  the  BOT 
surface  (Sg)  (Figure  4).  The  edge  currents  are  defined  as 


3*(P) 


tee 

J4 


'  tse 
J4 


p  e  S, 


p  e  S 


(amperes) 


(6) 
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BOT  edge 
patch 


Cap  edge 
patch 


P"' 


where 


-Triangle  function  for  edge 
transition  current,  Je 


Figure  4.  Detail  of  edge  region. 


j£e  -  ipfp(t)hp(P)  , 


J®e  -  ±  u  fz(t)hz(z),  z  »  +  L  , 
JC  z 


and  hp(p)  and  hz(z)  denote  half-triangle  functions  on  the  cap  and  BOT  edge 
regions,  respectively.  (The  t-components  of  the  currents  in  the  edge  region 
are  contained  in  the  expansion  set  for  Jct  and  58t.) 


4.3  Impedance  Expressions 


The  expansions  for  the  currents  on  the  BOT  surface,  the  caps,  the  wire 
segments,  the  junction,  and  the  edge  region,  given  in  Equations  (2)-(6),  are 
substituted  into  Equation  (1).  The  inner  product  of  the  integral  operators 
L(*)  in  Equation  (1)  is  formed  with  the  testing  functions  W  (■  1*)  via  the 


3 


The  right  column  vector  in  these  equations  represents  a  generalized  voltage 

vector,  and  the  Z  elements  are  the  familiar  impedances  defining  the  EM 

interactions  between  various  parts  of  the  body  in  Figure  2.  For  example,  the 

BOT-surface  interactions  are  given  by  Z8S  .  Similarly,  the  cap-cap,  the  BOT- 

mn 

cap,  the  wire-wire,  the  BOT-wire,  BOT- junction,  wire- junction,  and  the 

cc  SC  ww  sw  s  1 

junction- junction  interactions  are  defined  by  Zc  ,  Z  ,  Z  ,  Z  ,  Z  , 

m  mm 

Zw^,  and  Z^J,  respectively.  The  remaining  elements  refer  to  the  edge 
interactions. 

The  composition  of  the  system  matrix  in  Equation  (7)  is  determined  by  the 
particular  configuration  under  analysis.  For  example,  if  the  radiating  (or 
parasitic)  wire  element  in  Figure  2  is  not  attached  to  the  BOT  surface  (i.e., 
as  in  the  case  of  radiative  coupling),  the  junction-related  elements  can  be 
deleted  from  the  system  matrix,  Equation  (7).  Similarly,  if  there  are  no  off- 
surface  (wire)  antennas  on  the  BOT,  the  system  matrix  contains  only  elements 
of  Z",  Zcc,  Z®c,  and  the  edge-related  elements. 

Specific  analytical  expressions  for  the  various  matrices  in  Equation  (7) 
can  be  obtained  by  evaluating  the  inner  products  having  the  form 


Z  =  < 


w,  L(J)  >  =  JJ ds  JJ  ds'  jnk|ft*5  -  (7*w)(7'.j)j 


(8) 


where  all  superscripts  and  subscripts  on  W  and  J  are  omitted  for  simplicity 
and  S  and  S'  refer  to  surfaces  containing  the  field  and  source  points, 
respectively;  n  =  Ai/e  and  k  =  2n/X.  (For  the  choice  of  basis  functions  used 
in  this  discussion,  Equation  (8)  is  an  alternate  expression  to  the  one  given 
in  Equation  (la)  and  is  easier  to  program.)  Equation  (8)  is  evaluated  for 
each  of  the  regions  of  Figure  2.  Computer  implementation  forms  for  the 
impedance  elements  are  derived  in  Appendix  A. 


4.3.1  Junction-Independent  Impedance  Matrices 

First,  the  matrices  in  Equation  (7)  independent  of  the  junction  will  be 
considered.  The  derivation  of  the  elements  of  these  matrices  is  rather 
involved  and  will  be  relegated  to  Appendix  A.  Only  the  formal  definitions  and 
summary  of  results  are  given  here. 

s  s 

The  BOT  interaction  matrix  Z^  can  be  partitioned  into  submatrices 
corresponding  to  the  t  and  z  directed  current  components,  i.e., 


'  ss,tt | 

7ss,tz 

mn  | 

mn 

zS8,Zt| 

7ss,zz' 

_  mn  | 

mn 

(9) 


The  (i,j)-th  element  of  these  submatrices  is  obtained  from 


zss ,a0 
mn.ij 


L 

s 


(10) 


where  a,  0  are  combinations  of  t  and  z,  and  m  and  n  are  the  mode  numbers  of 

the  axial  expansion  (Equation  (2b));  wsf  and  are  the  i-th  testing  and  the 

mi  n  j 

j-th  current  basis  functions,  respectively,  on  the  BOT  surface. 

S  3 

Z  is  canonically  similar  to  the  system  matrix  arising  from  the  MM 
mn 

analysis  for  BORs  (Reference  5).  However,  in  the  present  case,  the  set  of 

expansion  functions  used  is  not  orthonormal  to  L  (•)  as  in  the  BOR  case,  and 

s 

3  S 

thus  there  is  no  modal  decoupling.  Hence  all  Z  matrices  are  present,  not 

mn 
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just  the  diagonal  ones.  Using  the  expansion  for  J  in  Equation  (2)  and  noting 
that  for  any  vector  A  in  the  BOT  coordinate  system, 


7*A 


3At 

3t 


3Az 

3z 


(11) 


s  s 

explicit  expressions  for  the  Z  elements  can  be  obtained  from  Equation 

mn 

(10).  These  matrix  elements  are  summarized  in  Appendix  A.l. 

The  cap  interaction  matric  Zcc  is  partitioned  into  t  and  p  (-  polar 
coordinate)  components  depicted  in  Figure  3,  i.e., 


Z 


cc 


zCC,tt  I  zCC,tp 

zcc,pt  J  zcc,pp 


(12) 


with  the  (k,  £)-th  element  defined  as 


cc  ,a£S 

k-X. 


+ca 
<Wk  , 


L 

c 


)>  , 


a,  S  *  t  or  p 


(13) 


where  W^01  and  are  the  k-th  testing  and  £-th  current  basis  functions  on  the 
corresponding  patches,  respectively.  For  maximum  generality,  the  patches  can 
be  trapezoidal.  For  any  vector  A  in  the  cap-coordinate  system. 


V*A 


A  3Ap  A^ 

p  3t  3r  r  » 


(14) 


where  r  =  rp  with  0  <  p  <  1  and  r  is  a  sector  ray,  i.e.,  the  distance  from  the 
origin  on  the  cap  to  a  sample  point  defining  the  generating  curve  of  the  BOT 
(Figure  3).  The  scaling  factor  r  may  be  different  in  each  patch.  Using 
Equations  (3)  and  (14),  the  expressions  for  Z^,a^  can  be  derived  as  shown  in 
Appendix  A. 2. 

s  c 

The  BOT-cap  interaction  matrix  Z  is  partitioned  as 

m 


zSC,tt  I  z8C,tp 

m  |  m 
zsc,pt  |  zsc,pp 


(15) 
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where  the  t  coordinate  on  the  caps  and  on  the  BOT  surface  is  the  same.  The 
(i,  4)-th  element  of  the  matrix  in  Equation  (15)  is  given  by 


-sc.ag 

m,i4 


■  <aj.  lc  <}“>»  . 


a  =  t  or  z;  8  *  t  or  p 


(16) 


where  WS<*  and  are  the  i-th  testing  function  on  the  BOT  and  the  4-th 
mi  4 

current  basis  function  on  the  cap,  respectively.  Explicit  expressions 

SC 

for  Z  are  derived  in  Appendix  A. 3.  The  cap-BOT  matrices  are  obtained  from 
„cs,a§  „sc,aB 

4*  n  .  *  u  .  n  • 

m,4j  -m,j4 

The  (k,  4)-th  wire-wire  interaction  matrix  elements  are  defined  as 


7ww 

k4 


<W*.  L(J^)> 


(17) 


where  is  the  k-th  wire  testing  function  and  is  the  4-th  wire  current 

ww 

basis  function.  The  explicit  form  for  Z^,  given  in  Table  1,  is  computed 
using  a  program  based  on  the  approach  in  Reference  11. 

The  BOT-wire  interaction  matrix  Zsw  is  also  partitioned  into  t  and  z 

m 

components,  i.e., 


zsw,t 

zsw  »  — 11 - 

m  Zsw,z 
m 


(18) 


with  the  (i,  4)-th  element  defined  as 


7sw,a 

n,i4 


<*s.  \  <*;»• 


a  ■  t  or  z 


(19) 


where  W8°  is  the  i-th  BOT  surface  testing  function  and  is  the  4-th  wire 
ni  4 

current  basis  function*  The  explicit  forms  for  the  Z  ) _  matrix  elements  are 

n,iJt 

derived  in  Appendix  A. 4.  The  wire-BOT  interaction  matrix  elements  are 

obtained  from  Zwa'a.  =  ZSW,<J..  An  analogous  development,  described  in  Appendix 
n,4j  — n,j4 

A. 5,  yields  the  cap-wire  interaction  matrix  elements,  Zcw. 
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TABLE  1.  JUNCTION  AND  WIRE  IMPEDANCE  ELEMENTS. 


4.3.2  Junction-Dependent  Impedance  Matrices 

The  derivation  of  junction-dependent  impedance  matrices  parallels  the 
development  in  Reference  3  for  junction  effects  on  BOR  surfaces.  For 
completeness,  a  synopsis  of  this  development  together  with  its  extension  to 
BOT  geometries  is  given  below.  The  junction  element  is  considered 
first.  Referring  to  Equation  (5),  the  junction  basis  function  consists  of  a 
disk  term  and  a  wire  attachment  term.  The  planar  annular  disk  patch 
approximates  the  BOT  surface  at  the  junction  point  (Figure  2c).  This 
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assumption  is  satisfactory  when  the  disk  radius  b  is  much  less  than  the  BOT 
radius  of  curvature  R  at  the  attachment  point  (i.e.,  b  «  R)  and  the  disk 
diameter  is  small  compared  to  a  wavelength.  In  practice,  the  latter  condition 
is  always  met  since  the  disk  diameter  is  taken  to  be  the  width  of  a  BOT 
surface  triangle  function  (i.e.,  usually  <  0.11  in  width).  Analytical 
predictions  using  this  approach  correlate  well  with  measured  data,  confirming 
the  efficacy  of  the  planar  disk  approximation.  Explicitly,  can  be  written 
as 


Zjj  «  <$-*,  L,(Jj)>  +  <wj,  L.(jj)>  +  <W^,  L.(Jj)>  +  <Wj,  L  (Jj )>  .  (20) 

a  J  a  djd  dja  ajd 

The  resulting  four  terms  describe  the  wire  (attachment)  segment,  the  disk,  the 
disk-wire,  and  the  wire-disk  interactions.  The  inner  products  in  Equation 
(20)  are  evaluated  using  Equation  (8)  and  the  fact  that 


V«3  = 


— i —  T 
2ita  a 


(h)  if  3  =  3J 

a 

H 


(21) 


-f  -f.  4 

and  similarly  for  v»w.  The  divergence  for  J  »  j-  is  expressed  in  terms  of  the 

d 

coordinate  system  of  the  planar  disk  and  implies  that  in  the  disk  region,  the 
radial  component  of  current  Jr  is  much  greater  than  the  azimuthal  component. 

In  evaluating  the  integrals  in  Equation  (8),  those  integrals  corresponding  to 
the  wire  attachment  segment  reduce  to  one-dimension.  Explicit  forms  for  all 
terms  of  Z^  in  Equation  (2)  are  summarized  in  Table  1. 

The  wire- junction  matrix  elements  are  defined  as 


ZV  ■  <*!•  vJi» +  <*?•  Lj<Jd)>  • 


(22) 


where  3^  is  the  i-th  wire  testing  function.  The  expressions  for  Z^  are  given 
in  Table  1.  The  junction-wire  matrix  Z^w  is  the  transpose  of  Zw^. 

The  BOT-junction  impedance  elements  are  given  by 


z‘jf  -  +  <c-  Lj(Id)>  • 


(23) 
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where  a  =  t  or  z,  m  is  the  mode  number,  and  $8“  is  the  testing  function 

Oil 

associated  with  the  i-th  BOT  axial  strip.  Explicit  forms  for  the  terms  in 

Equation  (23)  are  obtained  in  Appendix  A. 6.  The  junction-BOT  matrix  elements 

are  obtained  from  Z^8’a  *  Z8^’®  . 

n,i  -n,i 


4.3.3  Edge-Dependent  Impedance  Matrices 

The  edge  transition  region  is  subdivided  into  a  region  on  the  BOT  and  cap 
surfaces,  respectively.  (This  procedure  is  analogous  to  treatment  of  the 
junction  region.)  We  consider  first  the  interaction  between  the  k-th 
and  i-th  edge  regions,  which  in  terms  of  the  respective  testing  and  current 
expansion  functions  is  given  as 


Ci  “  <  «ke’  Le(JT>  >  +  <  «ke’  v5r>  > 


+  <  *r*  Leale)>  +  <  *k*.  Le(5J*)>  . 


(24) 


The  first  term  on  the  right  side  of  Equation  (24)  is  identical  to  the 
expression  for  z£^pp  (Equation  (A.16d)),  except  the  triangle  functions  are 
replaced  by  half  triangle  functions  denoted  as  h  functions.  The  remaining 
terms  of  Equation  (24)  are  evaluated  in  Appendix  A. 7. 

The  interaction  of  the  BOT  surface  with  the  edge  region  is  formally  given 
by 


7se,a 

m,U 


<*£.  >■.  <*T»  +  le(3r>>-  ‘  or  *> 


(25) 


where  the  first  term  on  the  right  side  is  Identical  to  Z8C'?P  where  the 

m,it 

triangle  functions  over  the  primed  variables  are  replaced  by  half-triangle 
functions,  i.e.,  the  g  functions  in  Equations  (A-20c)  and  (A-20d)  become  h 
functions.  The  explicit  expression  for  the  second  term  is  developed  in 
Appendix  A. 7. 

Similarly,  the  interaction  of  the  caps  with  the  edge  region  is  expressed 
as 


7ce,a 

Zki 


Le(}"):>  +  <5'“,  Le(J"» 


(a  ■  t  or  p)  . 


(26) 
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The  first  term  above  on  the  right  side  is  identical  to  when  the  g 

functions  subscripted  by  q  in  Equations  (A-16c)  and  A-16d)  are  replaced  by 
half  triangles.  The  second  term  is  derived  in  Appendix  A. 7.  All  edge 
interaction  matrices  retain  the  usual  symmetry  so  that,  for  example,  the  edge- 
BOT  interactions  are  given  by 

es,a  _se,a 
m,ii  -m,ii 

and  the  edge-cap  interactions  are  the  transpose  of  z£*,a. 

Finally,  the  wire  radiator-edge  interactions  are  defined  as 


*1- 


ce 


)>  +  ofj', 


se 


V3] 


)>, 


(27) 


where  the  first  term  on  the  right  side  is  identical  to  Z^j,p,  where  fj  is 
replaced  with  half  triangles  hj.  The  second  term  is  given  in  Appendix  A. 7. 
The  edge-wire  interactions  are  transposes  of  Equation  (27). 


4.4  General  Structure  and  Properties  of  the  Zp^.  Matrix 

The  Zg0T  matrix  in  Equation  (7)  is  composed  of  partitioned  submatrices 

g88,aB  (Equatjon  corresponding  to  the  orthonormal  vectors  spanning  the 

mn 

BOT  surface.  The  overall  structure  is  shown  in  Figure  5.  In  general,  this 
matrix  is  full  since  there  is  no  modal  decoupling.  (In  the  MM/BOR  analysis, 
only  the  matrices  Zmm  along  the  principal  diagonal  are  present.)  However, 
certain  symmetries  exist  for  and  that  reduce  the  fill  time  of  the  2jqT 
matrix.  Specifically, 


G  -  (G  ) 
mn  t  mn 

zsstt  „  {zsstt} 
mn  t  ran 
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t  mn 
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n  <  0  n  »  0  n  >  0 


where  t(Gmn)  denotes  the  transpose  of  G^.  In  the  implementation  of  the 
analysis,  only  the  lower  triangular  quadrant  of  the  ZgQ^.  matr^x  *-8  computed, 
i.e.,  the  partitioned  submatrices  are  filled  for  0  <  m  <  NMODE,  and 
-m  <  n  <  m,  where  NMODE  is  the  total  number  of  axial  modes  (including  m  ■  0) 
used  in  this  analysis.  The  matrix  fill  times  for  each  impedance  element  in 
MM/BOR  and  MM/BOT  are  about  the  same.  Furthermore,  detailed  examination  of 
the  Gmn(i,j)  function  shows  that  it  is  maximum  when  m  ■  n  and  i  ■  j,  i.e.,  the 
largest  values  occur  on  the  main  diagonal  of  the  ZBq<j.  matrix,  and  the  self- 


i 


terms  contribute  the  most.  These  properties  lead  to  a  diagonally  strong 
overall  matrix.  To  compute  the  inverse  of  ZgQ^,  the  entire  matrix  must 

be  filled,  which  can  be  accomplished  from  the  following  symmetry  relations: 


ntm 


zsstt  _  z88tt 
n,m  m,n 
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Z88tz  -  (Z88Zt) 

m,n  t  m,n 


sszz  sszz 
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These  relations  are  exploited  in  the  computer  algorithm  to  minimize  the  fill¬ 
time  and  aid  in  solution  of  the  matrix  equations  arising  in  the  formulation. 
If  the  BOT  surface  has  several  (physical)  planes  of  symmetry,  additional 
relations  can  be  established  within  each  of  the  partitioned  submatrices, 
permitting  further  computational  savings.  In  summary,  the  matrix  is 

composed  of  diagonally  strong  matrices,  possessing  certain  symmetries  and 
dominantly  intercoup led  only  for  neighboring  modes. 
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S.  FAR-FIELD  RADIATION  AND  SCATTERING  ANALYSIS 


Having  obtained  expressions  for  all  elements  in  the  overall  system  matrix 
(Equation  (7)),  the  unknown  current  vector  I  is  determined  when  the  applied 
excitation  (voltage  vector)  is  specified.  (In  a  scattering  analysis,  the 
incident  field  is  specified.)  Formally  the  radiated  field  in  8  (-  6  or  $) 
polarization,  found  from  the  currents  I  on  the  entire  configuration,  is^ 

Er  -ifi*  «'3kr  [**1  HI.  »0) 

where  the  row  transfer  matrix  Rp  consists  of  contributions  from  the  BOT,  the 
caps,  the  wire  (radiators),  the  junction,  and  the  edge  regions.  As  an 
example,  the  i-th  transfer  matrix  for  the  BOT  is  defined  as 


<JS“ 
'ni  ’ 


(31) 


where  the  superscript  a  denotes  t  or  z,  and  B  denotes  the  0  or  <j>  polarization 

g 

of  the  radiated  field  Ep,  given  in  spherical  coordinates  as. 


+6  +r  8in  9r  cos  ♦r  +  z  cos 

E^  -  ug  e  ,  (32) 


where  the  field  point  of  measurement  is  at  0^,  4>f,  and  p  is  the  distance  to  a 
point  on  the  BOT  surface,  measured  from  the  origin  (see  Figure  6).  Expressing 
the  inner  product  in  Equation  (31)  explicitly,  the  transfer  matrices  are  given 
as 


<*r\ 


J'dt  J  dz  (V“g)fJ(t)v“(z)e^  , 

i  -L 


(33) 


where  ■  p  sin  cos  +  z  cos  Referring  to  Figure  6,  the  unit  vectors 

in  Equation  (33)  for  the  BOT  geometry  are  given  by 

*  +  ♦ 
u_  -  u  cos  v  +  u  sin  v  ,  u 
t  x  y  z 
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i  th  strip 


and 


+r  ♦  ♦ 

u„*  u  cos  6  cos  *  +  u  cos  0  sin  *  -  u  sin  6 

0  x  rTry  r  r  z  r 

♦r  *  ,  + 

uA  *  -  u  sm  <t>  +  u  cos  4  . 

<j>  x  Yr  y  Yr 

Using  a  centroid  approximation  for  the  t-integration,  a  pulse  approximation 
for  the  triangle  functions  f“(*),  and  Equation  (2b),  explicit  expressions  for 
the  transfer  matrices  are  obtained  as  summarized  in  Table  2.  Analogously,  the 
transfer  matrices  for  the  cap,  wire,  junction,  and  edge  regions,  i.e.,  Rc,  Rw, 
rJ,  and  Re,  summarized  in  Table  2  can  be  obtained.  Using  these  R-matrices, 
the  far-field  power  radiation  patterns  and  the  scattering  cross  sections  can 
be  computed. 
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TABLE  2  TRANSFER  MATRIX  ELEMENTS 


BOT  Transfer  Matrix  Elements 


_s  ,t0 

R  .  »  cos 

ni 


0  DC  ^  At  T  C  cos(v  -  $  )Ar 
r  n  P  P  P  r'  p 

P-1 


r8;20  -  -  sine 
ni  r 


■:E 


At  T  2  Ar 
P  P  P 


P-1 


R8.t4>  „  jjt  V''  At  T  1  sin( v  -  4>  )Ar 
ni  n  P  P  P  r  p 

P-1 


R8*2*  -  0 
n,i 

where 


jk(x  cos  $  +  y  sin  $  )sin  0 
A  ■  e  P  E  P  r  r 

P 


DC  -  2L  sinc(C) 

n 


D2  -  2L  {sinc(£)  ~  (“l)n  sine  (K  cos  0r)} 


£  ■  nx  +  K  cos  0^,  K  ■  kL 
Cap  Transfer  Matrix  Elements 


C08  0 


16 

ZA  cos  ( v  -  ip  )  fC  gC  Ar 
P  P  r  P  P  P 

P-1 


TABLE  2  TRANSFER  MATRIX  ELEMENTS  (Continued) 
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rC j p 0  ^  co8  g  \  a  cos ( $  -  4>  )fP  gP  A 

\  r  ^  p  vyp  rr  p  p  | 

p-1 

16 

Rc,t(^  _  \  ''  A  S£R(V  _  $  )fc  gc  Ar 

k  p  p  p  p  p  p 

p-1 

16 


c ,  p<*> 


EA  sin($  -  $  )  fP  gP  Af 
p  vyp  yr'  P  6P  P 


p-1 


where 


ik  (x  cos  A  +  y  sin  A  )  sin  0  +  z  cos  0 

r  J  '  p  '  p  ~  V  IT  p  IT 

A  —  e 
P 


Mire  Transfer  Matrix  Elements 

4 

r”’0  *  /  {(Ax  cos  +  Ay  sin  $  )  cos  0  -  Az  sin  0  }  T  Ar 

l  /  j  1  p  yr  7p  r  r  p  rJ  p  p 

p-1 

4 

w,$  _  (-Ax  sin  i)>  +  Ay  cos  <p  )  T  Ar 

£  £,^4  P  r  p  r  p  p 

p=! 


where  Ar 
P 


is  evaluated  on  p-th  wire  segment 


Junction  Transfer  Matrix  Elements 


sin  a  cos  0  cos(v.  -  A  )  -  cos  a  sin 
q  x  d  Tr  q 
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TABLE  2  TRANSFER  MATRIX  ELEMENTS  (Concluded) 


RJ 2  8ln  °q  8in(Vd  '  V  Aq 


q-1 


where  Ar  is  evaluated  on  the  junction  disk  and  Rw,a  »  transfer  matrix  for 
q  a 

attachment  segment  Identical  to  RW,a  with  T  h 

P  P 

Edge  Transfer  Matrix  Elements 
8 


R*’6  “  ^  <Anf*(t)  h*(p)  A*  +  As  f*  (t)  h^(z)  A*’ 

i  L~i  P  P  P  P  PP  P  P 

P-1  < 


, *  *r, 

(ut*uB) 


where 


+  >r 

u  *u„  -  cos  0  sin(v  -  $  ) 
t  6  r  p  Yr' 


u  *u*  “  sin(v  -  <*■  ) 
t  $  P  r 


Ap  -  trapezoidal  edge  patch  on  cap 


As  -  edge  patch  on  BOT 
P 


and  A^  is  evaluated  at  the  center  of  Ap  on  the  cap; 

center  of  As  on  the  BOT. 

P 


Ar  Is  evaluated  at  the 
P 
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In  general  the  far-field  power  gain  in  B-polarization  is  given  by 


,  =  *LJL 

’6  4irP 


|R61 1 2 


(34) 


where 


Pq  =  Re{VZ_1V* } 


(34a) 


and 


I 


Z  lV 


(34b) 


The  current  vector  I  is  obtained  by  solving  the  system  of  linear  equations 
[Equation  (7)J,  subject  to  a  specified  voltage  excitation  V  determined  by  the 
location  and  electrical  characteristics  of  embedded  or  off-surface  wire 
antennas  on  the  body,  considered  in  more  detail  below. 


5. 1  Embedded  Antennas 

We  consider  first  the  case  of  embedded  (slot)  antennas  on  capped  BOT 
surfaces.  The  gain  is  calculated  from  Equation  (34),  where 


[RS,B  !  RC,B]  and  V  = 
n  I 


The  elements  of  the  R  matrices  are  given  in  Table  2;  I  is  computed  using  a 
system  matrix  Z  containing  only  the  impedances  of  the  BOT,  caps,  edges,  and 
their  respective  interactions.  It  is  assumed  that  there  are  no  antennas  on 
the  caps;  hence  Vc  »  0.  The  form  of  the  excitation  voltage  vector  V*  depends 


upon  the  type  and  Location  of  the  antenna  on  the  BOT.  As  an  example,  consider 
a  series  of  K  aperture  (slot)  antennas  embedded  in  the  BOT  surface  and 
centered  at  (t^,  z^),  k  =  l,  2  .  .  .  K.  Assuming  that  the  apertures  are 
rectangular,  then  if  the  aperture  at  the  i-th  axial  strip  subtends  one 
triangle  function  f^(t)  and  an  axial  width  of  (z^  -  Zq), 


In  the  above  expressions,  it  is  assumed  that  the  i-th  slot  is  excited 
uniformly  in  the  t  and/or  z  polarization  by  an  electric  field,  E^,  which  may 
be  a  complex  quantity.  In  the  program,  if  *  0,  then  (Ua)^  is  represented 
by  an  array  of  the  form 


(40) 


If  the  slots  are  excited  solely  in  the  t-polarization,  *  0.  For 

t  7 

simplicity,  and  can  be  set  equal  to  unity.  By  obvious  extension,  these 
results  can  be  generalized  for  a  BOT  containing  a  series  of  slots  at  a  given 
axial  strip,  as  in  the  case  of  two-dimensional  antenna  arrays. 

If  the  BOT  is  uncapped  (i.e.,  an  open  surface),  the  expression  in 
Equation  (34)  can  be  written  more  explicitly  as 


g 


B 


iSS 

mn 


(41) 


where 


P 

o 


V  Re  {Vs  Yss 
/  j  1  ra  mn 

ra,n 


(42) 
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ss 

The  matrix  Y  denotes  the  partitioned  submatrix  corresponding  to  the 
mn 

(m,n)-th  modes  in  the  inverted  Zgg-j,  matrix.  The  expression  for  the 

g 

matrices  remains  unchanged. 

5.2  Off-Surface  Radiators 

In  this  case,  the  gain  is  computed  from  Equation  (34)  where 


and 


the  closed  BOT  with  both  off-surface  wire  radiators  and  aperture 
and  VW  must  be  specified.  For  the  latter,  the  nonzero  element  of 
the  column  vector  corresponds  to  the  segment  of  the  wire  that  is  actively 
fed.  For  an  attached  wire  antenna,  also  is  specified.  For  a  parasitic 
radiator,  on  the  other  hand,  «  0.  Examples  using  this  formulation  for  wire 

radiators  are  given  in  Section  7.  In  all  cases  in  the  present  discussion,  we 
assume  that  Vc  and  V®  are  null  column  vectors,  i.e.,  no  active  radiators  are 
on  the  caps  and  edge  regions. 

5.3  Scattering  Fields 

In  general,  the  scattering  cross  section  in  terms  of  the  transfer 
matrices  in  Table  2  is  given  as 


n 

Vc 

Vw 

VJ 

„e 


If  we  excite 

antennas,  Vs 
n 


4xX2 


RP  Z  1Rq 


(43) 
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where  q  and  p  denote  the  polarization  of  the  incident  and  scattered  fields; 

(6  ,  *  )  and  (8  ,  4  )  are  the  incident  and  scattering  angles,  and  the  transfer 
i  i  a  s  s 

matrices  Rp  (g  -  p,q)  express  the  relations  between  the  current  on  the  BOT, 
m 

cap,  edge  regions,  and  an  observation  point  in  free  space.  Examples  of  cross 
sections  for  a  closed  cylinder,  represented  as  a  BOT  computed  using  this 
expression  are  given  in  Section  7.  In  general,  the  present  formulation  also 
can  treat  general  BOT  configurations  with  attached  wire  scatterers,  as  was 
done  for  BOR-wire  geometries  in  Reference  17. 


6.  NEAR-FIELD  AND  ANTENNA  COUPLING  ANALYSIS 


In  Che  preceding  sections,  Che  radiated  and  scattered  fields  were 
computed  at  field  points  sufficiently  distant  from  the  bodies  that  the 
respective  fields  were  planar  and  their  components  transverse.  Using  the 
previously  computed  currents  on  the  body,  the  near-fields  (electric  and 
magnetic)  at  points  <  \  distant  from  the  radiating  or  scattering  surfaces  can 
be  computed.  Generally,  all  six  field  components  must  be  determined  to 
specify  the  near-fields.  In  the  subsequent  sections,  the  present  formulation 
is  specialized  first  for  the  electric  near-fields,  followed  by  the 
corresponding  results  for  the  magnetic  fields.  Last,  antenna  coupling 
problems  will  be  addressed. 

6.1  Electric  Near-Field  Formulation 

In  general,  the  electric  field  at  a  free-space  point  r,  resulting  from  a 
surface  current  density  J  on  the  surface  of  a  body  in  Figure  2,  is  given  by 

E(r)  »  L  (JS)  +  L  (JC)  +  L  (j")  ,  (44) 

sew 

where  L(*)  denotes  the  integral  operators,  defined  earlier  in  Equations  (la) 
and  (lb),  over  the  currents  on  the  BOT,  cap,  and  wire  segments  of  the 
configuration  in  Figure  2.  The  contribution  from  the  junction  and  edge 
regions  is  assumed  negligible  and  Is  deleted. 

In  this  discussion,  the  current  is  induced  on  the  body  by  an  incident 
wave  (as  in  a  scattering  problem)  or  by  an  active  aperture  or  antenna  attached 
to,  or  radiatively  coupled  to,  the  BOT  surface.  It  is  convenient  to  reexpress 
the  operator  Lg(vfs)  in  Equation  (la)  for  this  analysis  as 


L  (J8) 

8 


< v- J® )  ( 


T  -  ?')  #jj 


ds' 


(45) 


using  the  relation 
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//< 


(V'»J)*(r  - 


r')  ds 


(V'.J)  (r  - 


r') 


ds' 


where  $^(r  -  r')  =  [(1/ R^)  +  j(k/R)}#(r  -  r')  and  $(r  -  r")  is  the  free-space 

Green's  function.  The  operator  for  L  (Jc)  is  identical  to  Equation  (45), 

c 

except  S  spans  the  trapezoidal  patches  on  the  caps.  Similarly,  L  (J*)  is  the 

w 

one-dimensional  analog  of  Equation  (45). 

First,  the  near-field  contributions  from  the  currents  on  the  BOT  surface 
alone  will  be  computed.  Substituting  the  expansion  for  J  (Equation  (2))  into 
Equation  (45)  yields 


£s(r) 


I8‘  J8‘ 
nj  nj 


*  Tsz 
u  I  « 
z  nj 


*(r'  -  r) 


3t 


?)*x(?'  -  h 


(46) 


where 


Using  the  pulse  approximations  for  the  triangle  expansion  functions. 
Equation  46  can  be  written  explicitly  as 


Eg(r)  “  jkn 


n,j  L 


nj 


At  (tc  gc 
£-j  p\  p  p 

P-1 


cos  v  +  —z  (x  -  x)  T 


T*  G*) 
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x) 


T2 

P 
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EAt  ( T 1  Gl  sin  v  +  (y  -  V)  TC  G1) 
P  V  P  P  p  k2  vyp  p  p/ 

p-1 


jmr  sz 

Vl  "J 
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E 

p-i 


y) 


Tz  G* 
P  P 


+  a 

z 

E 

n,j 

i8‘ 
n  j 

4 

E 

p-i 

-.&) 

TC  HC 
P  P 
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+  i  4 

nj 

4 

E 

p-i 

At  f] 

P  V 

z  z 

r  g 

p  p 

.  jnx  _z 

+  9  Tn 

k2L  P 

4\ 

(47) 


where 


L 


dz'  va(z')#(r  -  r)  ,  a  ■  t  or  z 

n  p 


L 

-■■I 


dz'  v  (z')  (z  -  z)#,(r  -  r)  , 

n  p  1  p 


(47a) 


(47b) 


and  G®  is  obtained  by  replacing  ♦  with  in  Equation  (47a).  The  G  and  H 
functions  are  evaluated  by  numerical  (Simpson)  integration  or  by  using  a 
Taylor  expansion  for  the  integrals. 

An  analogous  expression  is  obtained  for  the  near-fields  resulting  from 
the  currents  on  the  caps.  Replacing  5s  with  3  and  using  the  expansion  in 
Equation  (3)  yields 
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Ec(r)  *  jkn 


Vj-  IF- 


_ct  -ct  ct  .  _cp  ,cp  cp 
1l  f  A  8A  UpXA  f  A  gA 


']‘(?p  - 


b 


$  [*r  *?  •?  »;* +  iCte  *r  *r]  -  ?>v?P  -  f>)  •  <«> 


where  g^p  -  (r.)  *g»P  +  r  *gjp  an<l  S.  is  the  A-th  trapezoidal  patch  on  the 


’A  A  6A  '  A  “A  A 

caps.  The  surface  integral  is  evaluated  using  a  four-pulse  approximation  for 

the  triangle  functions  f  and  g.  Thus  for  the  A-th  patch. 


16 

ds'  ‘  H  \  ■ 


p=l 


where  Ap  is  the  area  intercepted  by  the  pulses.  The  total  field  attributable 


to  the  caps  is  obtained  by  summing  the  contributions  over  all  patches.  The 
resulting  field  can  be  decomposed  into  cartesian  components  if  desired. 


The  near-field  contribution  of  wire  radiators  (active  and/or  passive)  on 
or  near  the  BOT  can  be  obtained  by  using  the  expansions  for  the  wire  currents: 


E  (r) 
w 


jkn  S  /  d*  |Va  *(?A  -  ?)  +72  IT  (?A-?)*l| 
e  i-th  '  K 


(49) 


A  A-th  „ 
segment 


Using  a  centroid  approximation  for  the  integrals  and  a  four-pulse 


representation  for  j“ 

Xf 


I  T 

♦  V  ^  w  I  \  ^  -►  «► 

E  (r)  *  jkn  7  I,  (  u  /  Ax  T  *(r  -  r) 
w  J  >  A  )  x  P  P  P 


*  I  P-1 

4 


K  Y  AA  (x  -  x)T  •  (r„  -  r) 
k2  L-i  P  P  pip 

p-l 


+  h  1 


I  ♦1,1' 


(50) 


where  (x,  x  )  are  replaced  by  (y,  y  )  and  (z,  z  )  in  the  u  and  u 

r  r  r  y  Z 

components,  respectively.  If  the  wires  are  attached  to  the  BOT,  the  near¬ 
fields  include  the  junction  coupling  effects  since  i'*  is  computed  with  the 
entire  system  matrix,  including  junction  effects.  (The  direct  contribution  of 
the  junction  to  the  near-fields  is  negligible  and  was  deleted  in  Equation 
(44).) 

6.2  Magnetic  Near-Field  Formulation 

In  general,  the  magnetic  field  is  given  in  terms  of  the  vector  potential 
X  and  the  incident  field  *(.)  as 

fi(r)  -  frV)  +  VxX(r),  (51) 

where  r  is  an  arbitrary  test  point  where  the  field  is  to  be  sampled.  Letting 
«i(r)  -  0,  since  there  is  no  field  at  r  except  that  caused  by  currents  induced 
on  the  scattering  or  radiating  surface  S,  Equation  (51)  becomes 


(52) 


where  J  encompasses  the  contributions  from  the  BOT,  the  cap,  and  the  wire 
surfaces.  Again  the  direct  contributions  from  the  junction  are  omitted. 

The  near-field  contribution  from  the  BOT  currents  alone  is  computed 
first.  Noting  that  the  surface  current  density  on  the  BOT  can  be  decomposed 
into  t  and  z  components,  Equation  (52)  can  be  written  explicitly  as 


Using  the  modal  expansion  of  the  current  components  [Equation  (2)]  and 
evaluating  the  surface  integrals  in  the  same  manner  as  in  the  electric  near¬ 
field  analysis,  yields  the  following  expression  for  H  (r): 

s 


4 

+  “*  XI  Xnj  S  Atp((yP  '  y)  C0S  VP  '  (X 

n,j  L  P3*1 


] 


(54) 


where  the  y-component  is  identical  to  the  x-component,  with  the  substitutions 
sin  ♦  cos  and  (y  -  y)  +  (xp  -  x).  The  G  and  H  functions  are  defined  in 
Equations  (47a)  and  (47b). 

Similarly,  the  contribution  to  the  near-fields  from  the  cap  and  wire 
currents  are 
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Hc(r)  = 
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(55) 


where  the  y-component  is  identical  to  the  x-component,  with  sin  +  cos  and 
®lp  =  -  r).  Finally, 
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(56) 


where  in  the  y-coraponent,  (z„  -  z)Ay  ♦  (x  -  x)Az  and  (y„  -  y)Az  ♦ 

P  P  p  p  P  P 


(zp  -  z)Axp; 


in  the  z-component,  the  corresponding  terms  are  changed  to 


(yp  -  y)Axp  and  (xp  -  x)Ayp,  respectively. 


In  the  foregoing  discussion,  the  electric  and  magnetic  near-fields  were 
sampled  at  a  point.  An  alternate  approach  also  was  considered  in  which  the 
fields  were  sampled  and  averaged  over  a  rectangular  patch.  This  latter 
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formulation  was  a  generalization  of  that  given  In  Reference  31.  The 
predictions  from  the  patch  and  point  near-fleld  analyses  were  compared  for  a 
number  of  BOT  configurations.  Numerical  simulation  showed  that  In  regions 
where  the  electromagnetic  wave  departs  significantly  from  being  planar  (i.e., 
in  the  vicinity  of  apertures),  the  point-sampled  fields  were  more  accurate, 
and  thus  this  formulation  was  implemented  here. 


6.3  Coupling  Analysis 

The  concept  of  coupling  between  adjacent  antennas  can  be  expressed  in 
terms  of  the  modal  expansions  and  the  admittance  (or  impedance)  elements  of 
the  radiating  body.  The  definitions  and  results  obtained  are  a  generalization 
of  the  lumped  circuit  concepts  developed  in  Reference  36  (Chapters  13  and 
14).  In  general,  the  mutual  coupling,  defined  on  a  current  basis  (admittance) 
between  ports  i  and  j  for  a  multi-port  network,  is 


current  at  base  of  antenna  (port)  1 _ 

open  circuit  voltage  at  antenna  (port)  j 


(57) 


when  all  other  antennas  are  shorted  (i.e.,  *  0,  all  k  *  j).  This 

definition  is  specialized  below  for  coupling  between  wire  and  aperture 
antennas. 


Case  1.  Wire  Antennas  on  a  BOT  -  The  coupling  between  wire  antennas  on  a  BOT 

is  given  in  terms  of  the  junction  admittance  matrix  elements  of  the  respective 

wire  radiators.  Hence  if,  for  example,  a  BOT  has  two  monopole  wire  radiators, 

the  coupling  between  the  two  feed  points  is  given  by  Y^,,  i.e.,  the  junction 
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admittance  between  junctions  2  and  1.  Reciprocity  requires  that 

The  coupling  between  other  wire  radiators,  such  as  loops,  can  be  obtained  in 

similar  fashion. 


Case  2.  Aperture  Antennas  on  a  BOT  -  In  this  case,  we  must  consider  spatially 
averaged  coupling  coefficients.  An  example  will  illustrate  the  point. 

Consider  two  slot  antennas  (t-polarized) ,  each  subtending  single  triangle 
functions.  Let  the  slots  be  centered  at  (tA,  zA)  and  (tg,  zg)  and  have  an 
axial  half-length  of  LA  and  Lg,  respectively.  Equation  (57)  can  be  written  as 
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(58) 


<m“>-  <Ib> 

BA  ^  ’ 
A 


where  Che  spatially  averaged  currents  over  the  aperture  area  Sg  are 


‘S>  ■  ff  £  ^BA>mn  £<V 


ds  , 


SB  m,n 


(58a) 


and  similarly, 


•  fflji  £<ca>  ““ 


S .  m 
A 


(58b) 


Carrying  out  the  integrations  over  z, 


(y“)  sine 
BA  mn 


(=£) 


exp( jmirLg^L)  J  fftg)  dt 
tB 


(59a) 


=  XI  Vm  SlnC  (~ L~~ )  eXp(jmltLA/L)  /f(tA}  dt  ’ 


(59b) 


In  the  foregoing  expressions,  (Y„,)  is  the  element  of  the  admittance  matrix 
[inverse  of  the  system  matrix  in  Equation  (7)],  corresponding  to  the  axial  BOT 
strips  at  tA  and  tg,  and  the  m  and  n  modes.  By  analogous  considerations,  a 
similar  expression  can  be  obtained  for  <M^>,  associated  with  z-polarized 
slots. 
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7.  VALIDATION  AND  APPLICATION  OF  THE  GENERALIZED  MM/BOT  FORMULATION 


The  generalized  MM/BOT  formulation  was  usdd  to  predict  the  far  fields 
radiated  by  embedded  and  off-surface  wire  antennas  on  various  BOT 
configurations.  In  this  study,  both  active  and  passive  wire  radiators  were 
considered  on  open  and  closed  (capped)  BOTs,  and  the  near-fields  were  computed 
in  the  vicinity  of  such  antennas.  Finally,  the  coupling  between  antennas  was 
considered.  The  results  from  the  analysis  were  compared  with  available 
published  data  supplemented  with  independent  MDRL  experimental  measurements. 
Because  of  the  proven  reliability  and  accuracy  of  the  MM/BOR  codes  developed 
by  Mautz  and  Harrington  at  Syracuse  University,^  a  significant  part  of  the 
validation  effort  consisted  of  comparing  results  of  the  MM/BOR  and  MM/BOT 
formulations. 

In  subsequent  examples,  the  far-field  gain  calculations  are  given  in 
arbitrary  units.  When  these  calculations  are  compared  with  experimental  data, 
the  latter  are  used  as  the  normalization.  The  results  for  the  scattering 
cross  sections  and  the  near  fields  are  given  in  absolute,  unnormalized  units. 

7.1  Validation  of  the  Far-Field  Analysis 

7.1.1  Embedded  (Slot)  Antennas 

Because  of  the  availability  of  validated  MM/BOR  codes  that  treat 
apertures  in  surfaces  of  revolution, the  far-field  radiation  patterns  were 
computed  first  for  antennas  located  on  BOR  geometries,  such  as  capped  and  open 
cylinders,  modeled  as  a  special  case  of  a  BOT.  As  an  example.  Figure  7 
depicts  a  slot  antenna  of  2.06X  length,  subtending  an  angle  of  45°  and  located 
on  a  cylinder  with  a  radius  and  length  of  0.216X  and  2.76X,  respectively.  To 
facilitate  comparison  of  the  MM/BOT  and  MM/BOR  results,  the  z-axis  slot 
is  ^-polarized.  The  results  from  the  two  formulations,  depicted  in  Figure  8, 
show  excellent  agreement.  For  the  MM/BOR  calculations,  seven  circumferential 
modes  were  used;  for  the  MM/BOT  analysis,  four  axial  modes  were  employed.  The 
pitch  (vertical)  and  roll  plane  power  radiation  patterns  are  plotted  in  linear 
power,  normalized  to  the  MM/BOR  results. 
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7.1.2  Off-Surface  Antennas 


The  treatment  of  off-surface  antennas  was  examined  extensively  for 
monopoles  on  cylindrical  structures.  Initially,  the  simplest  case  was 
studied,  namely  monopoles  on  open  cylinders.  The  validation  baseline  was  the 
predictions  usirg  the  MDRL-developed  MM/BOR-wire  analysis  and  experimental 
data.^  The  pitch  and  yaw  plane  patterns  in  0  and  <j>  polarization  are  shown  in 
Figure  9  for  a  quarter-wave  monopole  mounted  at  the  mid-point  of  an  open 
cylinder  having  ka  =  1.357  and  a  length  of  2.76X.  The  MM/BOT  results  shown  in 
Figure  9  are  normalized  to  the  MM/BOR-wire  results.  Results  from  the  two 
analytical  techniques  are  in  close  agreement.  In  Figure  10  the  corresponding 
results  are  depicted  for  an  active  and  a  passive  (parasitic)  monopole  located 
on  the  same  cylinder.  Again  there  is  excellent  agreement  in  the  shape  and 
location  of  the  major  pattern  lobes  and  nulls. 

Next,  a  capped  cylinder  of  length  =  IX  and  ka  =  0.46  with  two  attached 
quarter-wave  elements  was  considered  (Figure  11).  The  element  at  the  cylinder 
mid-point  was  a  base-fed  monopole;  the  other  attached  element  was  unfed 
(passive).  The  MM/BOT  predictions,  compared  with  the  results  in  Reference  3, 
are  summarized  in  Figure  11.  The  radiation  patterns  were  computed  with  and 
without  the  effect  of  the  junction  (current)  interaction  at  the  base  of  the 
unfed  element.  The  predicted  and  measured  results  in  Figure  11  clearly 
indicate  the  importance  of  including  the  junction  interactions.  In  this  case, 
the  unfed  stub  acts  as  a  reflecting  element  dramatically  altering  the 
directivity  of  the  fed  element.  (The  result  without  the  junction  interactions 
models  a  nonattached,  radiatively  coupled  element,  and  the  pattern  can  be 
obtained  using  the  analysis  in  References  11  and  33.)  The  patterns  are 
plotted  in  power  using  a  linear  (not  decibel)  scale  and  normalized  to  unity. 
Nine  circumferential  modes  were  used  in  the  MM/BOR-wire  analysis  and  four 
axial  modes  for  the  MM/BOT  results. 

An  application  of  the  present  formulation  to  wing-mounted  antennas  is 
depicted  in  Figure  12a  where  the  yaw  plane  patterns  for  a  monopole  mounted  on 
the  trailing  edge  of  a  wing  section  are  depicted.  The  section  is 
2.76X  long.  The  wing  coordinates  are  given  in  Sample  Problem  3,  Vol.  II  of 
this  report.  In  Figure  12b  the  effect  of  an  asymmetric  active  element  on  the 
radiation  pattern  is  shown.  Junctions  were  used  at  the  antenna  attachment 
points. 
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9  -  polarization 
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- MM/BOR 


6  -  polarization  <)>  ■  polarization 

Yaw  plane  (dB)  Yawp  lane  (dB) 


Figure  9.  Comparison  of  MM/BOR  and  MM/BOT  computed  power  radiation  patterns: 
cylinder-mounted  monopole. 
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6  ■  polarization 
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<t>  -  polarization 
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Figure  10.  Comparison  of  MM/BOR  and  MM/BOT  computed  power  radiation  patterns: 
cylinder-mounted  opposing  monopoles  (active  and  passive). 
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0  -  polarization 
Pitch  plate  (dB) 


Figure  11,  Comparison  of  MM/BOR  and  MM/BOT  computed  power  radiation 

patterns  with  experiment:  cylinder-mounted  active  and  passive  monopoles. 
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1/8  X  monopole  at  mid-plane 


8  ■  polarization 
Frequency  =  96.15  MHz 
Yaw  plane  (dB) 


Figure  12.  Power  radiation  patterns  for  wing-mounted  monopoles. 


7.2  Validation  of  the  Scattering  Analysis 


Examples  validating  the  scattering  analysis  In  Section  5.3  are  considered 
below.  In  Figure  13  the  blstatic  scattering  cross  section  for  an  open 
cylinder  of  radius  0.216A  and  length  2.76A  is  given.  The  cylinder  is 
illuminated  broadside  (6  *  90°)  with  a  TE  wave.  The  absolute  cross  sections 
predicted  by  the  MM/BOR  and  MM/BOT  formulations  shown  in  Figure  13  are  in 
close  agreement.  The  scattering  cross  sections  for  the  foregoing  cylinder 
with  caps  is  depicted  in  Figure  14  and  compared  with  experimental  data 
obtained  by  Oshiro  et  al.  '  The  inclusion  of  the  edge  region  terms  is 
critical  for  obtaining  satisfactory  results,  particularly  in  the  horizontal 
polarization. 

The  applicability  of  the  present  analysis  to  degenerate  BOT  surfaces  is 
demonstrated  in  Figures  15  and  16,  where  the  monostatic  scattering  cross 
sections  for  a  flat  plate  (2A  on  a  side)  and  the  bistatic  cross  section  for  a 
thin  disk  of  0.8A  diameter  are  depicted.  For  the  former  case,  the  MM/BOT 


2.76  A 


Figure  13.  Comparison  of  MM/BOR  and  MM/BOT  computer  bistatic  cross  section  for  open  cylinder: 
broadside  illumination,  TE  polarization. 
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calculations  using  4  axial  modes  and  16  triangle  functions  are  compared  and 
are  in  good  agreement  with  the  experimental  results  of  Ross.^  The  minor 
discrepancies  are  due  to  the  fact  that  the  experimental  data  were  measured  for 
a  relatively  thick  plate.  In  Figure  16,  the  bistatic  scattering  cross  section 
for  a  disk  is  shown,  comparing  the  MM/BOT  and  MM/BOR^  formulations.  The 
agreement  is  excellent  in  spite  of  the  trapezoidal  representation  of  the  disk 
in  the  MM/BOT  calculations.  The  disk  was  subdivided  in  the  MM/BOT  analysis 
into  24  patches. 
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°  Capped  MM/BOT 
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Horizontal  polarization 


Figure  14.  Measured  and  predicted  monos  tat  ic  cron  section  for  capped  cylinder. 
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Figure  15.  Monostatic  scattering  cross  section  for  a  square  plate. 


7.3  Validation  of  Near-Field  Analysis 

Validation  of  the  near-field  formulation  for  the  electric  and  magnetic 
fields  is  shown  in  the  subsequent  examples.  In  this  phase  of  the  investiga¬ 
tion,  the  fields  were  computed  by  the  MM/BOT  technique  for  points  at  distances 
«  1  from  the  radiating  or  scattering  surfaces.  Unfortunately,  there  is  a 
paucity  of  measured  near-field  data  suitable  for  comparison  with  the  present 
analysis.  As  one  benchmark,  the  exact  boundary  value  solution  for  an  infinite 
right-circular  cylinder  fed  with  a  ^-polarized  slit  was  used  (Figure  17).-^ 

The  slit  subtended  an  angle  of  deg  and  was  excited  with  a  uniform  field. 

For  fields  sampled  near  the  cylinder,  the  predictions  of  the  MM/BOT  analysis 
for  a  finite  cylinder  can  be  compared  with  the  exact  (classical)  solution 
obtained  for  the  infinite  cylinder.  A  comparison  of  these  solutions  for  a 
cylinder  with  ka  =  1.35  and  a  slit  of  45°  is  shown  in  Figure  18.  For  the 
MM/BOT  analysis,  the  cylinder  length  was  chosen  to  be  2.761  and  the  slit 
length  was  21.  The  near  fields  are  sampled  along  a  line  bisecting  the 
aperture,  resulting  in  »  0.  The  closest  field  sampling  was  at  1.2a,  which 
corresponds  to  0.2a  (0.041)  from  the  plane  of  the  aperture.  The  calculations 
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—  —  —  ox*/X21  MM/BOR 

_ _____  „x8  n2  I  (Mautz  and  Harrington) 

J  2  MM/BOT 

o  oxa  f)r  (using  trapezoidal  patches) 


Figure  16.  Comparison  of  MM/BOR  and  MM/BOT  calculated  bistatic  cross  section 
for  a  thin  conducting  disk;  a*®  »  E-plane  pattern,  o'*'*-'  -  H-plane  pattern. 

using  four  and  seven  modes  produced  practically  the  same  results.  Also  shown 
in  Figure  18  are  the  fields  obtained  from  a  patch  near-field  formulation  in 
which  the  sampled  fields  are  averaged  over  a  flat  strip.  The  large 
discrepancy  of  the  patch  results  from  the  exact  solution  at  points  near  the 
body  occurs  because  in  this  region  the  EM  wave  departs  significantly  from 
being  planar.  At  distances  ~  10A,  the  patch  and  point  formulations 
coalesce.  (As  expected,  at  these  distances,  the  classical  solution  for  the 
infinite  cylinder  and  the  BOT  results  for  the  finite  cylinder  diverge.) 

In  Figure  19,  the  near-field  components  sampled  at  a  radial  line  at 
^  *  45°  to  the  aperture  center  are  given.  Again,  the  results  of  the  exact  and 
the  MM/BOT  solutions  are  in  excellent  agreement.  The  corresponding  results 
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Figure  17.  Classical  solution  for  slit  cylinder  (0-excited  slit). 


for  field  points  sampled  at  0  =  90°  are  shown  in  Figure  20.  At  sampling 

distances  »  X  from  the  BOT,  E  decreases  and  the  EM  wave  front  tends  to 

P 

become  planar.  In  all  these  calculations,  4  modes  were  used  and  17  points 
defined  the  circumference  of  the  cylinder.  As  a  final  example,  the  near-field 
was  computed  for  the  wing  configuration  in  Figure  12b.  The  result  is  given  in 
Figure  21,  showing  the  nonplanar  nature  of  these  fields.  Four  axial  modes 
were  used  in  this  analysis. 
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Field  intensity 


Figure  19.  Computation  of  near  fields  for  slit 


0-fed  slit  cylinder  (ka  =  1 .35) 
o  Classical  analysis 
cylinder  length  =  “> 
antenna  aperture  •  infinite  45°  slit 

a  MM/BOT :  point  near-field  analysis 
IMP  =  17 

NMODES  =  4  and  7 
cylinder  length  *  2.76X 
antenna  aperture  =  2  A  x  45° 


at  0  *  45°  Wit  angle  •  45°). 


1.2a  2a  10a  5X  10X 
Sampling  points 


0- fed  slit  cylinder  (ka  =  1 .35) 

o  Classical  analysis 
cylinder  length  =  °° 
antenna  aperture  -  infinite  45°  slit 

*MM/BOT :  point  near -field  analysis 
NP  =  17 
NMODES  =  4 
cylinder  length  =  2.76X 
antenna  aperture  =  2X  x  45° 


Figure  20.  Computation  of  naar  fields  for  slit  cylinder  at  <t> ■  90°  (slit  angle  *  45°). 


8.  COMPUTER  IMPLEMENTATION 


The  MM/BOT  formulation  described  in  the  preceding  sections  was  imple¬ 
mented  with  the  A-STAR  computer  algorithm  described  in  detail  in  Volume  II  of 
this  report.  The  overall  structure  of  the  program  flow  compares  with  that  of 
the  MM/BOR  codes  and  the  original  MM/BOT  code  (Reference  2).  The 
computational  complexity  is  approximately  equivalent  for  both  MM/BOT  and 
MM/BOR  for  a  given  size  body.  The  matrix  fill-times  are  comparable.  The 
major  differences  lie  in  the  fact  that  the  modes  in  the  present  analysis  do 
not  decouple  as  in  the  MM/BOR,  although,  in  general,  the  resulting  network 
matrices  remain  diagonally  strong  and  have  certain  symmetries.  The  off- 
diagonal  submatrices  for  m  *  n  used  in  the  MM/BOT  formulation  can  sometimes  be 
deleted  from  the  computation  without  excessive  error  penalty. 

Adequate  computational  accuracy  is  achieved  when  the  BOT  surface  and  the 
wire  radiators  are  segmented  into  strips  <  0.15A  in  width  or  length.  The 
number  of  axial  modes  chosen  is  dependent  upon  the  spatial  accuracy  desired 
for  the  surface  currents.  While  the  examples  shown  in  the  validation  section 
involved  mostly  right-circular  cylinders,  the  present  formalism  is  capable  of 
treating  any  asymmetric  BOT,  such  as  wing  sections  as  is  discussed  in  Volume 
II. 
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APPENDIX  A:  DERIVATION  OF  THE  IMPEDANCE  MATRICES 


The  derivation  of  the  matrix  elements,  contained  In  Equation  (7)  and 
discussed  In  Section  4,  Is  detailed  below. 


cc 

A.l  Impedance  Matrix  for  the  BOT  Surface,  Z 

— i - mn 

Using  Equation  (8),  the  elements  of  the  partitioned  submatrices 

of  ZSS  ate  computed  from 
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Using  the  expansion  for  Js  (Equation  (2))  and  the  expression  for  the 

divergence  in  Equation  (II),  and  noting  that  WS^  =■  (JS*)*,  it  can  be  shown 
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where 
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v  and  V*  are  the  angles  between  the  t-curve  and  the  x-axis  at  points  (x,  y) 
and  (x't  y')  on  the  BOT,  respectively.  Similar  expressions  can  be  obtained 


for  the  other  components  of  Z 


ss 

mn* 


The  triangle  function  f^(.)  and  its 


derivative  fj(*)  in  Equation  (A-2)  are  approximated  by  four  pulses,  denoted  as 
TP  and  ^p»  resPectively,  with  p  =  1,  2,  3,  4,  where 


P  v  v  i)  and  *P  =  l1*  _1*  -1!  * 


Using  the  functional  forms  for  v^(*)  [Equation  (2b)]  and  centroid 
approximating  the  t  and  t'  integrals,  yields 
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and  similarly, 
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where  the  function  Gmn(p,q)  is  the  Green's  function  integrated  over  ds  and 
ds',  where  these  integrations  can  be  carried  out  quasi-analytically,  as 
discussed  below,  and 
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In  the  above  expressions,  the  convention,  adopted  in  Reference  5,  is  used  such 
that  the  indices  (i,j)  correspond  to  the  indices  (p,q)>  respectively.  The 
corresponding  triangle  functions  for  these  indexes  are  depicted  in  Figure 
2b.  The  term  is  the  angle  between  the  t-curve  and  the  x-axis  at  point  (xp, 
yp)  on  the  BOT. 

Evaluation  of  Functions  -  The  evaluation  of  Gmn  encompasses  two  surface 
integrals,  written  explicitly  as 


g  (p,q) 

mn 


exp[jit(-  mz  +  nz")/L]  #  , 


(A-7) 


where  the  t  and  t'  integrations  are  carried  out  over  the  p  and  q-th  subseg¬ 
ments  (strips)  associated  with  the  i-th  and  j-th  axial  strips  on  the  BOT  (see 
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Figure  2b).  Letting  £  =  — = -  ,  the  two  integrations  over  z  and  z'  reduce  to 
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R  - 


(x  -  x")2  +  (y  -  y')^  ,  and  K  ■  kL. 


Evaluating  the  integrals  over  t  and  t'  by  centroid  approximation  over  the 

intervals  At  and  At  yields 
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pq 


where  Rp^  is  obtained  by  evaluating  p  at  p  ,  i.e.  , 


(A-l  1 ) 
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pq 


+  <y. 


An  alternate  evaluation  of  Equation  (A-8)  can  be  carried  out  by  using  a 
Taylor-expan8ion  about  tn  and  approximating  the  resulting  integral  by  a  series 

r 

as  is  done  in  Reference  5.  The  present  procedure  is  computationally  more 
streamlined  than  the  former. 

For  the  self-terms  (i.e.,  p  “  0),  Equation  (A-7)  reduces  to  the  form, 


G 

mn 


At 

jK  -jJ  +  (1  +  jKO 


The  integral  over  £  can  be  evaluated  numerically  via  Gausalan  quadrature  or  a 
Simpson  integration  routine  after  an  integration  hy  parts  la  used  to  soften 
the  In  £  singularity  at  the  lower  limit.  Thus, 
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A. 2  Impedance  Matrix  for  the  Caps,  2C 

Specializing  Equation  (8)  for  the  caps  yields 
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sk 


[a,  8  =  t  or  p  (=  polar  coordinate)] 


( A- 1 4 ) 


where  the  surface  integrals  are  over  the  k  and  £-th  patches.  Using  Equations 
(3)  and  (14),  we  obtain 
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(The  independent  variables  t  and  p  are  suppressed  in  f  and  g  for  notational 
simplicity.)  The  integrals  in  Equation  (A-15)  can  be  evaluated  using  a  four- 
pulse  representation  for  the  triangle  functions  fk  and  g^,  denoted  by  fp  and 
gp,  yielding 
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where  Gp^  denotes  the  Integrated  free-space  Green's  function  over  the 
trapezoidal  areas  subtended  by  the  p-th  and  q-th  pulses  and  corresponding  to 
the  k-th  and  A-th  patches,  respectively.  Similarly, 
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For  the  nonself  term  case,  k  *  A, 
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where  A^  and  A^  are  the  trapezoidal  areas  subtended  by  p-th  and  q-th  pulses, 
and  4  is  the  Green's  function  evaluated  at  the  center  of  each  area.  For  the 

pq 

self-terms,  k  -  A  or  p  *  q,  Gp^  is  evaluated  by  carrying  out  the  integration 
over  a  given  trapezoidal  patch  shown  in  Figure  3,  yielding 
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where  r  “  y  (r^  +  The  lflSt  integration  is  done  numerically. 


A. 3  Impedance  Matrix  for  the  BOT-cap  Interactions,  Z 
For  this  case,  Equation  (8)  becomes 
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(A-18) 


where  and  S'  denote  the  i-th  axial  strip  on  the  BOT  and  4-th  patch  on  the 

caps,  respectively.  Using  Equations  (2)  and  (3)  for  the  current  expansions 

and  Equations  (11)  and  (14)  for  the  divergence  terms  and  noting 

that  w8“  »  (J8®)*,  Equation  (A-18)  can  be  written  explicitly  as 
mi  mi 
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(The  independent  variables  in  the  f,  g,  and  v  functions  are  implied.)  As  in 
Sections  A.l  and  A. 2,  using  pulse  approximations  for  the  f  and  g  functions  on 
both  the  BOT  strip  and  the  cap  patch,  yields 
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where 


(A-21) 


and 
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G  (p.q)  is  evaluated  by  numerical  integration, 
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A. 4  Impedance  Matrix  for  the  BOT-Wire  Interactions.  Z 


In  this  case  Equation  (8)  becomes 
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where  the  integrals  are  over  the  i-th  BOT  strip  and  the  i-th  wire  segment,  the 
latter  being  a  line  integral.  Substituting  the  appropriate  testing  and 
expansion  functions  into  Equation  (A-22)  yields 


<C"l)n  ■  Jkn  ffd‘  /dh'  jv3«  *7  *;«>*>'>{  \ 


(z)  *  , 
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(A-23) 

where  u^  is  a  unit  vector  along  the  i-th  wire  segment.  Using  pulse 
approximations  for  the  triangle  functions  on  the  BOT  strip  and  wire  segments. 
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( A-24a) 


and  likewise, 


( A-24b) 


where  Gm(p,r)  is  defined  in  Equation  (A-21)  where  the  index  r  is  associated 
with  the  Jt-th  wire  segment  and  Aq  ♦  A^. 

A. 5  Impedance  Matrix  for  the  Cap-Wire  Interactions,  Zcw 
The  elements  of  Zcw  are  given  by: 


sr-3«  -  ?  <-sr>  1 


,w 


♦  ,  (a  *  t  or  p) 


(A-25) 


where  the  integrals  are  over  the  k-th  cap  patch  and  the  Jt-th  wire  segment. 
Substituting  the  testing  and  expansion  functions  into  Equations  (A-25)  and 
using  a  pulse  representation  for  the  triangle  functions  on  the  caps  and  the 
wire  yields 
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where  gP  =  (r  )  *  gp  +  r  ^g£,  Ap  is  the  area  of  the  trapezoidal  patch  spanned 
by  the  pulse  functions  on  the  caps,  and  *  is  the  free-space  Green's 
function. 


A. 6  BOT-Junctton  Impedance  Matrix,  Z*J 

The  BOT-junction  impedance  elements  consist  of  a  wire  segment  part  (first 
term  of  Equation  (23))  and  a  disk  part.  The  former  are  identical  to  the  BOT- 
wlre  interaction  matrices  (Equation  (A-24))  except  that  the  integration  £s 
over  the  attachment  segment  spanned  by  a  half-triangle  function  (i.e.,  in 

ra  1 

Equation  (A-24)).  Using  Equation  (5b)  in  Equation  (8),  the  disk  part 
of  2s-1  for  the  t  and  z  components  can  be  obtained  as  follows: 
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where  is  over  the  i-th  axial  strip  of  the  BOT  and  s'  is  over  the  disk 
region  (Figure  2c),  i.e., 
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Using  Equation  (2b)  and  a  pulse  representation  for  f,  the  foregoing  integrals 
can  be  approximated  as 
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where 
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(A-29) 
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and  $  is  the  free-space  Green's  function  with 
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(xp  -  xq)2  +  (yp  -  yq)2  +  (z  -  zq)2 


The  integration  over  at  is  treated  by  a  centroid  approximation  of  the  integral 
at  four  equispaced  angles  on  the  disk,  mid-distant  between  the  outer  and  inner 
rims. 


A. 7  Edge-Dependent  Impedance  Matrices 

in  this  section,  the  edge-edge,  BOT-edge ,  cap-edge,  and  wire-edge 
impedance  elements  are  derived. 

The  edge-edge  elements  are  formally  given  by  Equation  (24).  The  first 
term  on  the  right  side  of  Equation  (24)  is  identical  with  zj^,PP  in  Equation 
(A-16d),  except  that  the  f  functions  are  represented  by  half  triangle 
functions  h,  and  are  denoted  as  ^^,PP»  The  second  term  incorporates  the 
edge-edge  interactions  on  the  BOT,  i.e.. 
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where  for  notational  brevity,  the  independent  variables  are  suppressed  and  the 
z-superscripts  on  f  and  h  are  implied  in  this  discussion.  Using  a  pulse 
representation  for  f  and  h, 
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where  AS  =  At  Az  and  4>  is  the  free-space  Green's  function  evaluated  at  the 

p  p  p  pq 

center  of  each  pulse.  The  edge  interaction  between  the  BOT  and  the  caps  is 
similarly  given  by 


<Ke’  \  <jr» ■ 


A  fZ  hZ  fP  hP  $ 

q  p  p  q  q  pq 


(A- 30b) 


where 


hp  =  (7  )-1  hp  +  r"1  hp  . 

q  q  q  q  q 


Following  the  convention  of  the  previous  sections  of  the  appendix,  the 

indices  (p,  q)  refer  to  (k,  i) ,  corresponding  to  the  unprimed  and  primed 

coordinates,  respectively.  The  edge  interactions  between  the  cap  and  the  BOT 

are  the  transpose  of  Equation  (A-30b).  The  total  edge-edge  interaction 

between  the  k-th  and  i-th  segments  is  a  sum  of  the  foregoing  expressions 

and  Z^*PP. 
k£ 

Following  the  approach  above,  the  interactions  of  the  BOT  surface  with 
the  edge  region,  formally  given  by  Equation  (25),  can  be  written  as 
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where 
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and  it  (z)  is  the  free-space  Green's  function  with  R  (z)  =  j (x  -  x  )  + 
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Similarly,  the  interactions  of  the  cap  surface  with  the  edge  region, 
given  in  Equation  (26),  can  be  written  explicitly  as 
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(A-32b. 


where  is  defined  as  hP  before. 

P  P 

Finally,  the  wire  ( radiator)-edge  interactions  between  the  i-th  radiator 
and  the  j—th  edge  segments,  in  Equation  (27),  can  be  written  in  terms  of  pulse 
representations  of  the  triangle  functions  as 
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